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■rf , ABSTRACT 

o 

' Bose-Einstein correlations and invariant momentum distributions are analyzed 

for longitudinally expanding finite systems, like jets in elementary particle colli- 
sions or systems created in high energy heavy ion reactions. 
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1. Introduction 



p Bose-Einstein correlations are in general not measuring the whole geometrical 

1^ ■ sizes of big and expanding finite systemJlH'l since the expansion may result in strong 

correlations between space-time and momentum space variables not only in the lon- 
gitudinal, but in the transverse and temporal directions, tool. 
X : Where have all the geometrical sizes gone? One can sho Jli, that they are dis- 

\ guised in the invariant momentum distribution of the bosons in case they cancel from 

the radius parameters of the Bose-Einstein correlation function (BECF). 

We shall briefiy review herewith the results presented in refs.Bl'i§ which shall be 
appended with an application to jet size determination. 

2. Wigner Function Formalism 

The two-particle inclusive correlation function is defined and approximately ex- 
pressed in the Wigner function formalism as 

C(Ak:K) = —^5 — ^ — -~1H = ^. (1) 

(n)' iVi(pi)iVi(p2) I 3(0, K) |2 ^ ^ 

In the above line, the Wigner-function formalismlli'B is utilized assuming fully chaotic 
(thermalized) particle emission. The covariant Wigner-transform of the source den- 



sity matrix, S{x,p) is a quantum-mechanical analogue of the classical probability 
that a boson is produced at a given x'^ = (t, r) = (t, r^., r^, r^) with = {E,p) = 
{E,px,Py,Pz)- The auxiliary quantity S{Ak,K) = J d'^xS{x, K) exp{iAk ■ x) appears 
in the definition of the BECF, with Ak = pi — P2 and K = {pi +p2)/2. The single- 
and two-particle inclusive momentum distributions (IMD-s) are given by 

, , E da ~ , , , , EiEo da , , 

N^{p) = —— = S{Ak = {),p), and iV2(pi,p2)= , , , (2) 

otot dp atot dpi dp2 

where atot is the total inelastic cross- sect ion. Note that in this work we utilize the 
following normalization of the emission functionB: / -^d'^xS{x,p) = (n). 

3. Effects from Large Halo of Long-Lived Resonances 

If the bosons originate from a core which is surrounded by a halo of long-lived 
resonances, the IMD and the BECF can be calculated in a straightforward manner. 
The detailed description is given in ref.0, here we review only the basic idea. 

If the emission function can be approximately divided into two parts, representing 
the core and the halo, S{x; K) = Sc{x; K) + Sh{x\ K) and if the halo is characterized 
by large length-scales so that ShiQmin] K) « SciQmin] K) at a finite experimental 
resolution of Qmin > 10 MeV, then the IMD and the BECF reads as 

iVi(p) = iVi,,(p) + iVi,,(p), (3) 
C{Ak-K) = l + K^S^!^^, (4) 



where A^^i j(p) stands for the IMD of the halo or core for i = h,c and 



A* = \^{K = p) 



iVi(p) 



H 2 



(5) 



Thus within the core/halo picture the phenomenological A* parameter can be obtained 
in a natural manner at a given finite resolution of the relative momentum. This 
parameter has been introduced to the literature by Deutschmann long time agoi. 
In the core/halo picture, the effective or measured intercept parameter A(p) can he 
interpreted as the momentum dependent square of the ratio of the IMD of the core to 
the IMD of all particles emitted. 

4. General Considerations and Results 



We are considering jets in elementary particle reactions or high energy heavy ion 
reactions, which correspond to systems undergoing an approximately boost-invariant 
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Figure 1: Emission of particles with a given momentum is centered around and 77^ 
for systems undergoing boost-invariant longitudinal expansion, as indicated by the 
shaded area. 

longitudinal expansion. For fully boost-invariant longitudinal expansions, the emis- 
sion function may depend on boost-invariant variables only. These are defined as r = 
r^ = Q.hH {t + z)/{t-z)i mt = ^E^-pl y = 0.5 ln[(£; + / - p,)] 

and Tt = + r^. For finite systems, the emission function may depend on i] — 
too, where yo stands for the mid-rapidity. Approximate boost-invariance is recovered 
in the | rj — yo \<< Ay region, where the width of the rapidity distribution is denoted 
by Ay. In terms of these variables the emission function can be rewritten as 

Sc{x;K)d'^x = Sc,^{T,7],r^,ry)dTTodr]dr^dry. (6) 

The subscript * indicates that the functional form of the source function is changed, 
and it stands for a dependence on K and Un also. 

In the standard HBT coordinate systemEl, the mean and the relative momenta are 
K = {Ko,Kout,0,KL) and Ak = {Qo,Qout,Qsidf^^L)- Note that the side component 
of the mean momentum vanishes by definitionBS. Since the particles are on mass- 
shell, we have = K ■ Ak = KqQq — K^Ql — KoutQout- Introducing (3l = Ki/Kq and 
jSout = Kout/Ko, the energy difference Qq can thus be expressed as 

Qo = PlQl + PoutQout- (7) 

If the emission function has a such a structure that it is concentrated in a narrow 
region around {Ts,ris) in the (t, 77) plane, then one can evaluate the BECF in terms 



of variables r and 1] by utilizing tlie expansion 

Ak-X = Qot- QoutTx - QsideTy - Q z ^ (8) 

Qt^ Qout^x Q side^y - Qn'^siv - Vs)- (9) 
Tfie coefficients of tlie r and the rs(?7 — r^^) are new variables given by 

Qt = Qo cosh[i]s] - Ql sinh[?7s] = {PtQout + PlQl) cosh[r7s] - Ql sinh[?7s], (10) 
Qr} = Ql cosh[r7s] - Qo sinh[?7^] = Ql cos1i[?7^] - {(3tQout + (^lQl) sinh[r7s]. (U) 

In terms of these new variables the BECF reads as 

C{/\k- K)C^I+ ^ 1 + K{K) \S'^^*^Q^Q^^QouuQs^de) \\ 

|S(0,ir)|2 I S,,,(0,0,0,0) |2 

At this level, the shape of the BECF can be rather complicated, it may have non- 
Gaussian, non-factorizable structure. Gaussian approximation to eq. (|T2|) may break 
down as discussed in more detail in the Appendix of ref.i. 

5. Mixing Angle for HBT 

In the experimental analysis, one of the most frequentlyS but not exclusively0 
applied parameterization of the BECF is some version of a Gaussian approximation. 
The out-longitudinal cross-term of BECF has also been discovered in this context 
recently^. In order to identify how this term may come about, let us assume that 

Sc,*{r,r],r:^,ry) = H^{t) G^{r]) I^{r^,ry). (13) 

In Gaussian approximation one also assumes that 

H,{r) oc exp(-(r-r,)V(2Ar2)), (14) 
GM « exp(-(r/-r/,)V(2Ar/2)), (15) 
I*{rx,ry) oc e-XY>{-{{rx- rx,sf + {ry- ry^sf)/{2Rl)). (16) 

The corresponding BECF is given by a diagonal form as 

C{M- K) = l + X^ exp(-Q2 Ar^ - Qlr^Ar^l - Q]RI). (17) 

This diagonal form shall be transformed to an off-diagonal one if one introduces the 
kinematic relations between the variables QriQ-q and the variables QoutiQi- In the 
HBT coordinate systemll one finds 

C{Ak-K) = i + KeM-RldeQlde-RLtQLt-RlQl-^Rlut,LQoutQL)m 

Rlide = Rl^ (19) 



R-lut - -R* + ^R-lut^ (20) 
6Rl, = P^icosh'[ns]AT: + smh'[r^,]T^Ar^l), (21) 
Rl = ipL smh[r]s] - cos\i[r]s]fTlAril + {[3l cos\i[r]s\ - sm\i[r]s]f At"", , (22) 
Rlut,L = (A cosh [t],] {13 l cosh [r],] - sinh [q,])) Ar^ + 

(A sinh[r/,] (/?i sinh[r7,] - cosh[r/,]))r2 Ar^.l (23) 

Note that the effective temporal duration, Ar* and the effective longitudinal size, 
TgArj^ appear in a mixed form in the BECF parameters SRl^^, R\ and R^out li 
their mixing is controlled by the value of the parameter ru- These results simplify a 
lol in the LCMS, the Longitudinally Co-Moving Systemll^l, where /3l = 0: 

5i?L = PK^os\?[ris]ATl + sm\i\,^sysAril), (24) 
Rl = cosh2[r7,]r2Ar/2 + sinh2[r/,]Ar2, (25) 
Rlut^L = -Asinh[r/,]cosh[ry,](Ar2 + r2Ar/2). (26) 

Let us define the Longitudinal Saddle-Point System (LSPS) to be the frame where 
Vsi'iTT't) = 0. In LSPS one finds that 

5Rl, = P^AtI (27) 
Rl = TlAr,l + f3lATl (28) 
= PthArl (29) 



Introducing Qo = PtQout + PlQl and Qt = \/Qlut + Qlide ^^e BECF can be rewritten 
in LSPS as 

C{Ak-K) = l + KeM-^rlQl-rlAr,lQl-RlQl). (30) 

Thus the out-long cross-term can be diagonalized in the LSPS frameBi. The cross 
term should be small in LCMS if rj^^^^ « 1 i.e. if | y — yo |<< Ayi. Since the size 
of the cross-term is controlled by the value of t]s in any given frame, it follows that 
r]s is the cross-term generating hyperbolic mixing anglS for cylindrically symmetric, 
longitudinally expanding finite systems which satisfy the factorization of eq. ([T3|) . 

6. Decaying Lund Jets 

An ultra-relativistic jet corresponds to a boost-invariant longitudinally expanding 
finite system, with strong correlation among longitudinal space-time and momentum- 
space variables. Finite correlation lengths are created by the random fluctuation of 
the break-up points of the hadronic string. Local left-right symmetry of the jet 
fragmentation prescribes the following proper-time distributiontJ: 

H{T)dT = — -6^+"(/€r)2"+i exp(-6(fi:r)2)/trfr (31) 

r(i a) 



where k ^ 1 GeV/fm, a = 0.3 and b = 0.58 GeV ^ are the default parameters of 
JETSET7.4 0. In Gaussian approximation to the BECF one obtains 

1 r(3/2 + a) , . If r2(3/2 + a)\ , , 

''s= r) = ^= ^ / , and Ar^ = il + a - \/ / • 32 
V6k2 r(l + a) bn \ T {1 + a) J 

The BECF for prompt pions may become measurable as discussed in ref. 0. If the 
prompt pions cannot be separated, the above expressions need to be corrected for the 
resonance decay effects which may increase both the mean and the variance of the 
distribution. The transverse radius parameter for direct pions is Rside = -R* ~ Rstring 
and the correlation length between space-time rapidity and momentum space rapidity 
is to be evaluated numerically. This can be performed bvutilizing the Lund mapping 
of momentum space to space-time along the lines of refs.ll3llj. When resonance decays 
were switched off, one obtainedEl Arj^ ~ 0.5 using JETSET6.3 default parameters of 
a = 0.5 and b = 0.9 GeV~^. This yields Rl ~ ToAr/,, ^ 0.6 fm and Ar* ^ 0.5 fm/c for 
prompt pions. When resonance decays are switched on, the width of the (r/ — y) 
distribution (for all pions) increased to Arj^ ^ 1 unit. Should the measured BECF 
for direct pions be fitted with the expression 

C{Ak, ir) ^ 1 + A. expi-QlAr^ - QIrI - Q?i?L»,), (33) 

the calculated numbers could be contrasted to data and some information on the 
transverse size of the string may also become available. 

It should be clear that even a successful description of the BECF for particles 
from a decaying string, like the famous Andersson-Hofmann mo del0 does not directly 
reveal the total longitudinal size of the string because the BECF is dominated by the 
correlation length between space-time and momentum space rapidityi. 

In ref.i one of us argued that the total longitudinal sizes of expanding finite 
systems may become measurable with the help of a combined analysis of IMD and 
BECF measurements. According to ref.i, the space-time rapidity distribution of the 
boson source can be measured as the asymptotic large transverse mass limit of the 
invariant momentum distribution. 

These results may be considered as first steps into a new direction of jet size 
determination with the help of combined IMD and BECF measurements. 



7. Geometrical vs. Thermal Length Scales for Heavy Ions 

For high energy heavy ion reactions, we model the emission function of the core 
with an emission function described in detail in ref.i. This corresponds to a Boltzmann 
approximation to the local momentum distribution of a longitudinally expanding 
finite system which expands into the transverse directions with a transverse flow, 
which is non-relativistic at the maximum of particle emission. The decrease of the 
temperature distribution T{x) in the transverse direction is controlled by a parameter 



a, the strength of the transverse flow is controlled by a parameter b. A parameter 
d controls the strength of the change of the local temperature during the course of 
particle emissioni. If all these parameters vanish, a = b = d = 0, one recovers the 
case of longitudinally expanding finite systems with T{x) = Tq with no transverse 
flow, as discussed in ref.i, ifa = d = Oj^b the model of refS is obtained. 

The parameters of the correlation function are related by eqs. ( p!8| - p9D to the 
parameters R^,, and TsArj^ which in turn are given by 

^2 = ^ + ^cosh[r7,], (34) 

Here the geometrical sizes are given by Rq, the transverse size, At], the width of 
space-time rapidity distribution and At, the duration around the mean emission 
time Tg = Tq. The hyperbolic mixing angle 77^ at midrapidity yo i, where also the 
out-long cross-termll3 vanishes. The thermal length-scales (subscript t) are given by 

RT = -r-l^T7^ ^Vt = tt, Att = -^ — . (37) 



The transverse mass of the pair is denoted by Mf = y-R'o ~ -^1- 

These analytic expressions indicate that the BECF views only a part of the space- 
time volume of the expanding systems, which implies that even a complete measure- 
ment of the parameters of the BECF as a function of the mean momentum K may 
not be sufficient to determine uniquely the underlying phase-space distribution. 

It is timely to emphasize at this point that the parameters of the Bose-Einstein cor- 
relation function coincide with the (rapidity and transverse mass dependent) lengths 
of homogeneit'^ in the source, which can be identified with that region in coordi- 
nate space where particles with a given momentum are emitted from. The lengths 
of homogeneity for thermal models can be obtained from basically two type of scales 
referred to as 'thermal' and 'geometrical' scales. 

The thermal scales originate from the factor exp(— p • u{x)/T{x)), where u{x) 
is the four-velocity field. Figure 2 and 3 illustrate how the temperature changes 
in the transverse or temporal directions induce transverse mass dependent thermal 
radius or thermal duration parameters. This is to be contrasted to the 'geometrical' 
scales, which originate from the exp(/i(a;)/T(x)) factor which controls the density 
distributioni. Here fi{x) stands for the chemical potential. 
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Figure 2: The temperature gradients in the transverse direction create a transverse 
mass dependent effective thermal radius parameter. In this illustration, the effective 
size of the region where particles with a given rrit are emitted from is decreasing with 
increasing values of rrit. Note that the transverse flow gradients may result in a similar 
effect^'^, not indicated on this illustration. 




Figure 3: Temporal changes of the local temperature create a transverse mass 
dependent effective thermal duration, which is decreasing with increasing values of 
rrit in this illustration. 



8. Dynamically Generated Vanishing Life-Time Parameter 



As a consequence of the possibility for temporal changes of the local temperature, 
we find that the effective duration of the particle emission Ar^, becomes transverse 
mass dependent and for sufficiently large values of the transverse mass this param- 
eter may become vanishingly small. The reason for this new effect is rather simple: 
Particles with a higher transverse mass are effectively emitted in a time interval when 
the local temperature (boosted by the transverse ffow) is higher than the considered 
value for m^. If the local temperature changes during the course of particle emission, 
the effective emission time for high transverse mass particles shall be smaller than 
the effective emission time of particles with lower transverse mass values. 

For a more detailed analysis of the model the interested reader is referred to ref.iJii, 
where it is pointed out that under certain conditions the parameters of the Bose- 
Einstein correlation function may obey an Mj-scaling: Rgide — Rout — Rl oc I/a/M^. 
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